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^ ; Abstract 

\ The current accelerated universe could be produced by modified gravitational dynamics as it can be seen in 



particular in its Palatini formulation. We analyze here a specific non-linear gravity-scalar system in the first 



constant problem exists also in this formalism, in close parallel with the standard metric formulation. 



^ ■ order Palatini formalism which leads to a FRW cosmology different from the purely metric one. It is shown 
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' that the emerging FRW cosmology may lead either to an effective quintessence phase (cosmic speed-up) or 

00 : 

to an effective phantom phase. Moreover, the already known gravity assisted dark energy dominance occurs 
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^ , also in the first order formalism. Finally, it is shown that a dynamical theory able to resolve the cosmological 
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I. INTRODUCTION 



The dark energy problem (for astrophysical data indicating dark energy existence and related 
universe acceleration see and refs. therein) is considered to be a challenge for modern Cosmol- 
ogy. There are various approaches to its construction. In general, one can classify most efforts as: 
a) modified gravity (for instance, the so-called f{R) theory); or: b) introducing some exotic matter 
with negative pressure. The truth, as usual, may lie in between. Indeed, since the two directions 
above are completely phenomenological one can go further and try to modify gravity via a non- 
standard (non-minimal) coupling with matter. This may in fact provide the effective description 
for a dark energy universe. The importance of such models has already been pointed out in P| 
where a possible origin for dark energy dominance was suggested. In P, Q a suitable non-linear 
matter-gravity theory was used to construct the dynamical approach to the resolution of the cos- 
mological constant problem. (Such coupling may be used for cosmological purposes as in ;5]). Here, 
the gravity-matter coupling part is assumed to be a f{R) theory-like Langrangian, non-minimally 
coupled with a scalar field Lagrangian. It is interesting that in such models the standard conformal 
transformation to scalar-tensor theory is over-complicated (if possible at all) . This implies that in 
the purely metric formalism only one possible frame really exists. A first aim of ours is to see how 
much the picture changes if one allows a first-order a la Palatini approach to the gravitational field. 



In the recent study of modified gravity with a non linear 1/R term ja] (for its one-loop 
quantization see [3]) in relation with cosmic acceleration it was further drawn the important lesson 
that the metric formalism is not equivalent to the Palatini (first-order) formalism. (Actually, this 
fact was known in general after the early proof given in |3l]). Indeed, while modified gravity of only 



special form may be viable in the metric version |9| (due to a possible instability 10]) its Palatini 



version may be instead viable (see 1^, ll^l ^-nd references therein). Of course, the cosmological 
dynamics of both versions are in principle different. Notice, however, that modified gravity may 
be always rewritten as a (well-studied) scalar-tensor theory. 

It is clear that a generic non-linear matter-gravity theory may lead also to non-equivalent 
gravitational physics, depending on the choice of either metric or Palatini description. 
In 0| it has been shown how dark energy may become dominant over standard matter at small 
curvatures in the metric formulation. This is directly caused by the choice of the Lagrangian. 
Non-linear gravity-scalar field systems and modified gravity may nevertheless provide an effective 
phantom or quintessence dark energy with the possibility of explaining the present speed-up of 
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the universe (and without the need of introducing negative kinetic energy fields) 
It has been proven in that a dynamical approach to the cosmological constant based on a 
suitable non-linear gravity-scalar field system with a scalar field potential having a minimum at 
a generic, but negative value provides a possible solution to the cosmological constant problem. 
This is based on the presence of a non-standard kinetic term in the scalar field Lagrangian, 
which diverges at zero curvature. In this particular model the cosmological constant freezes at a 
nearly zero value, while field equations are stable under radiative corrections and they provide the 
expected asymptotic limit. 



The same class of non-linear matter-gravity theories are here studied in the Palatini formalism 
(where the purely gravitational part is suitably generalized) and applied to alternative gravity. This 
formalism, considering the metric and the connection as independent fields, produces second order 
field equations which are not affected by instability problems and seem to be in a fairly acceptable 
accordance with solar system experiments The Palatini version leads to a completely new 

formulation when compared with the standard metric formulation and it provides cosmological 
models which are able to explain also the present time cosmic speed-up. 

It is also shown that the Palatini formalism in the case of non-linear gravity-scalar systems gives an 
acceptable realization of the dark energy dominance and the resolution of the cosmological constant 
problem. Thus, even if the metric and the Palatini formulations apparently lead to different 
gravitational physics, one can see that the same cosmological phenomena may qualitatively occur 
in both frameworks, which therefore deserve to be carefully compared. 



II. FIRST-ORDER NON-LINEAR GRAVITY-MATTER SYSTEMS: 
THE FIRST ORDER FORMALISM 

We consider a 4-dimensional Lorentzian manifold M, with a metric g^^ and an independent 
connection T^^. The starting non-linear gravity-matter Lagrangian is: 

L = ^{F{R)+ f{R)Ld) + KL,nat{'^) (1) 

where f{R) and F{R) are some analytic functions of the scalar field R = g'^^ RapiT) and is a 
scalar field Lagrangian; here we set Ld = —\g^^d^(f) dycl) + V {(j)) = —\ \\ V0 p +¥{({)). In prin- 
ciple, Ld could be any other matter Lagrangian (describing spinors, vectors, etc). Moreover, one 
can think of more general cases when f{R) (or even a function of Ricci squared terms .14]) couples 
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to the Riemann (Ricci) tensor contracted with some derivatives of matter and multiphed by 
a similar matter Lagrangian. It is important to remark that in all cases (even if is chosen to 
be the Lagrangian of a scalar theory!) the general class of theories described by phenomenological 
Lagrangians of the form cannot be conformally transformed into an equivalent scalar-tensor 
theory, with the only exception of the theories already treated in 12:] and :13?] (i.e. all the theories 
in which = 1, i.e. f{R) theories). The term Lmat represents a matter Lagrangian, functionally 
depending on arbitrary matter fields ^ together with their first derivatives, equipped with a grav- 
itational coupling constant k = SttG, which can be supposed to be k = SvrG = 1 in natural units. 
In the first order Palatini formalism the Ricci-like scalar is R = R{g, T) = q°'^ Rnfi(T) where R^u{r) 
is the Ricci tensor of any independent torsionless connection F; see p, \\^ for details. Equations 
of motion are respectively: 

1. Field equations for the metric field g: 

F'{R)R^, - \F{R)g^, + /' {R) LaRf,, - \f (R) Lag^, + / (R) V = r^. (2) 

where T^i, = — d,y(j) is stress-energy tensor of the scalar field and r^i, = TJJIf^ is 
the matter stress-energy tensor generated by taking the variational derivative of Lmat{^)- 
Taking the trace of equation with respect to the metric g'^'^ we obtain the scalar-valued 
master equation: 

F'{R)R - 2F{R) + (/' {R) R-2f (R)) La + f{R)T = T (3) 

where r = g^^T^y and T = g^'^T^u = Ld — V {(p) = \\ V0 |p. The above equation can be 
re- written, more conveniently, as: 

2F{R) - F'{R)R + r={f'{R)R-f (R)) - f {R) V (0) 

Hence, the scalar- valued equation @ controls solutions of equation We shall refer to 
this scalar-valued equation as the structural equation of spacetime. For any real solution 
R = F{t,T) of © it follows that f{R) = f{F{T,T)) and f'{R) = f{F{T,T)) can be 
seen as functions of the traces r and T of the stress-energy tensors of matter and of the 
scalar field (p. For convenience we shall use the abuse of notation f{T,T) = f{F{T,T)) and 
/'(r,T)=/'(F(r,r)). 



Thus, in the same way as in 
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the generalized Einstein equations may be introduced: 



Rfiu (T) — g^aPy 
4 



(4) 



where one defines the operator PlJ" = — ^^^^Ty + \tu . The scalars b and c, depending 
on are here defined as: 

b = b{R) = F'{R) + f (R) La 

(5) 

c = ciR) = i {FiR) + / (R) La) = ^^^^^^ 



(see also Section llll|) . 

2. Field equations for the connection field T: varying the Lagrangian with respect to the 
connection one gets: 

Vl[^g^^''{F\R) + f{R)La)]=0 (6) 

which states that T can be chosen to be the Levi-Civita connection of the metric h^u = 
{F' (R) + /' (R) Ld) g^u = bgf_iu, providing the spacetime manifold with the typical bi- metric 
structure ensuing from the first-order Palatini formalism This implies that the gener- 
alized Einstein equation @ can be immediately rewritten as: 

R,.u (h) = g^aP^ (7) 

where R^i,{h) is now the Ricci tensor of the new metric /i, being T = Tic{h)- 

3. Field equations for the scalar (dilaton-like) field (j): performing the variation of the La- 
grangian with respect to (f) we obtain: 

du {^f{R)g^''d^^) = -^f{R)V' (0) (8) 



4. Field equations for the matter fields ^ : they can be obtained from the Lagrangian by 

mat ( 



the usual prescription ^^"^^1^^^ = 0, depending on the particular form of L„iat{'^)- 



Before proceeding further let us make two important notes. Firstly, in order to get Eq.(6) one 
has to additionally assume that L^at is functionally independent of an arbitaray connection F; 
however it may contain metric covariant derivatives of some matter fields denoted here by ^ . 
This means that the matter stress-energy tensor T^"* = r^^, depends solely on the metric g and 
fields ^ together with their partial derivatives. Thus, usually this tensor is covariantly conserved, 
in a sense of metric covariant derivative, i.e. V^r^jy = 0. Nevertheless, we do not assume a priori 
that the (metric) covariant derivative of the left-hand side of Eq. (2) identically vanishes (in fact 
it does not in general). However this does not provide any inconsistency in our model. When we 



solve the field equations of motion with a source - here the right-hand side of Eq. (2) - having it 
a vanishing covariant diveregence thus also the left-hand side of Eq. (2) taken on this particular 
solution (i.e. "on shell") will have automatically vanishing covariant derivative as well. 

Secondly, we are now close to discuss energy momentum conservation in modified gravity models. 
As it is well-know that any generally covariant field theory admit the covariant conservation of 
energy momentum which can be expressed in the form of generalized Bianchi identitiy. This 
general fact has been presented in many books and articles and it is widely recognized as (second) 
Noether theorem (see e.g. a recent book and ref.s quoted therein). More in particular, for 
the case of modified (including scalar-tensor) theories of gravity (also within a Palatini formalism) 
the issue has been nicely discussed in a recent paper , where explicit formulae for the Bianchi 
identity can be also found. 



III. FRW COSMOLOGY FROM GENERALIZED EINSTEIN EQUATIONS 



Owing to the astonishing cosmological experimental results recently obtained in 1] and to 
the renewed interest in alternative theories of Gravity for cosmological applications, it is quite 
interesting to study the cosmological implications of the formalism above. Let 5 be a Friedmann- 
Robertson- Walker (FRW) metric: 

1 



g 



-dr + a\t) 



(9) 



where a(t) is a scale factor and K is the space curvature (K = 0, 1, — 1). The stress-energy tensor 
of a matter (ideal fluid) is 



pmat 



= {P+ p)Uf,Uy + pg^y 



where p is the pressure, p is the density of matter and is a co-moving fluid vector, which in a 
co-moving frame (n^ = (1,0,0,0)) becomes simply: 



P 












o\ 






(10) 



pa^ity"^ 
\ pa'^{t)r'^sm^{e) ) 

The standard relations between the pressure p, the matter density p and the expansion factor a{t) 
are assumed 



p = wp 



p = rja 



(11) 



where particular values of the parameter w E {—1,0, |} will correspond to the vacuum, dust or 
radiation dominated universes. Exotic matter (which is up to now under investigation as a possible 
model for dark energy) admits instead values of w < —1; see e.g. jlj|. These expressions (|11|) follow 
from the metric covariant conservation law of the energy-momentum 'VgT^u = (as dicussed above) 
and consequently the continuity equation should hold: 

p + 3H{p + p) = (12) 

where H = - is the Hubble constant. 

a 

We recall here that the field equations for the metric field g can now be rewritten in terms of the 
new metric h, as already done in 0. Relying on the methods already developed in it is simple 
to see that the generalized FRW equations following from Q are 

where P/ and Pq are the components of Plf = — ^^y^T^ + \tI} . We recall here that the 

parameters b and c are defined ©• Owing to the cosmological principle and to the consequent 

homogeneity and isotropy of spacetirne, we argue that T^^,, t^^ and depend just on the time 

parameter t. As already noticed in |l2 . \\\ for much simpler cases, we get from the structural 

equation (see the above discussion) that both b and c can be expressed just as functions of time. 

, ■ 2 

In our case, we have moreover that Ld= 2't' + ^ {4')^ so that: 

pi _ C f{R) rpl , i_i _ c , P 

— b b ^1 ^ b^i — b ^ b 

(14) 

pO _ c f{R)rp0.in_c f{R)j^ p 
/ — 6 ~~b~^0 ^ b^O — b ~2b~^ b 

The modified FRW equation can be thus written as: 

and the non-minimally coupled scalar-field equation reduces to: 

j^{^f{R)4) = ^f{R)V'{<l^) (16) 

From now on we assume, for simplicity and following the headlines of P|, that V {(j)) = 0. This 
implies that the field equations can be written as ^f{R)(f) = const or equivalently gf{R)'^Lii = 
= const. It follows that, on-shell with respect to the above field equations for the scalar field 
<j), the non-minimal coupling term in the Lagrangian takes the form 

/(i?)Lrf = (3^f{R)-'a~^ (17) 




which differs from the results of |j| obtained in the metric formahsm in the particular case F{R) = 
R, which is in fact degenerate (see below). 



A. A particular case: F{R) = R 

We start here considering the degenerate case F{R) = R. This case is worth to be considered 
as it is the simplest case of non-linear gravity-scalar system reproducing General Relativity in 
the limit f{R) = 0; results can moreover be simply compared with results already obtained in 
P] in the metric formalism. We stress that to obtain an exactly solvable class of models we will 
hereafter consider in larger detail the particular cases / {R) = aR^ (a 7^ 0; n 7^ 1) for matter- 
free universe with vanishing spatial curvature {K = 0), since these cases are in fact analytically 
solvable. Nothing would prevent a qualitative study of more general cases. We will specify these 
conditions step by step in the following, when explicitly necessary. 
In the case F{R) = R the master equation gives: 

R + T = P^ f{R)-^ (/' {R)R-f (R)) (18) 

which simply becomes an (analytic) equation for R: 

R+{3w- 1) ,a-3(i+-) = (19) 

In the case of polynomial (or more generally analytic, if solvable) functions f{R), solutions R of 
the above equation are, in general, rational (analytic) functions of the variable a. Consequently, 
as already remarked before, R can be simply expressed as a function of time. For any solution 
R = R{a(t)) one can further calculate: 

6 = 1 + /' (i?) Lrf = 1 + (R) f{R)-^a-^ 

(20) 

c = \{R + f{R)Ld) = \{R + (3'' f{R)-'a-^) 
The r.h.s. of the generalized FRW equation becomes: 

ya 2b) a2 66/(i?)a6 ^36 66a3(i+'") ^ ' 

and it is consequently a rational function of a. 

Let us now consider the specific example where / {R) = aR^ (a 7^ 0; n 7^ 1). It is supposed that 
there is no matter content in these models. We thus obtain, from the structural equation: 

R = aR^ (n - 1) (22) 



so that we have aR^L^ = f (R) = ^^^-R- Apart from the trivial solution R = (which implies 
6=1, c = and + ^ = 0) there exists also a nontrivial solution corresponding to the case 
1 = aii"^^ (n - 1) Ld. In this non-trivial case the time-depending parameters are found to be 
(n 7^ i to ensure 6 7^ 0): 



b = l + f'{R)Ld = l + ^ 



2n-l 
n-1 



(23) 



c=liR + aR-Ld) = \R(l + ^^)= R^ 
The conformal parameter b is constant and the modified Friedmann equations are 

a) ^^~6 b ^36~ 6b ^ 6{2n-l) 6 \{n-l)b^ 2^^) ^ ^ 
which in turn can be re-written in the simplified form: 

(^V K l„n + l 

As we already stressed before, the coefficient b is now constant; then R = g^'^R^u (bg) = R (g) is 
the true Ricci scalar of the FRW metric g, which is known to be equal 

2' 



R{g) = 6 



a a 
-+ - 
a \a 



6H^{1- q) (26) 



where q is the deceleration parameter. Hence the two metrics g and h are equivalent, apart from 
a constant conformal factor. The generalized Friedmann equation (specified now to the exactly 
solvable case K = 0) takes the form: 

which immediately gives the value of the deceleration parameter 

^ = ^ = -l + A (28) 

without having to solve explicitly an ordinary differential equation for the scale factor a{t). 
It is however possible to obtain the explicit dependence of the scale factor on time. We recall 
that it has been assumed from the beginning (see (|16l) ') that V{(p) = and consequently the field 
equation for the scalar field is: 

|(V5i?» = (29) 

which in turn implies that y^i?"^ = const or equivalently that on-shell R^^L^ = f3'^a~^. 
From the structural equation (|2U() we consequently obtain, in the case under discussion, that 



r 91 

R = [(n — 1) a/5 J a "+i . Therefore, again under the exphcit further assumption i^T = 0, we 
have: 



n + 1 
6(2n- 1) 



[(n - 1) a0^] 



L(n, Q, /3)a^' 



3 

' n + l 



(30) 



where we have defined for simphcity the constant factor L(n, a, /3) 



71+1 



[(n- l)a/32] 2("+i). 



6(2n-l) 

Solving this simple differential equation we get the explicit expression of the scale factor as a 
function of time: 



a{t) 



3L{n, a, (3) 
n + l 



(31) 



where ts is some integration constant. More precisely, Friedmann equation in the simplest form 
(or equivalently a ~ a^~'^), where 7 = M^^ifLCtil^ implies that the deceleration parame- 



ter is g = 7 — 1. Moreover, in the case 7 > 0, choosing the solution corresponding to an expanding 
universe one has a ~ t^, i.e. a forever expanding universe with Big Bang at the origin tg = 0. 
Conversely in the case 7 < 0, changing the arrow of time in order to avoid shrinking solutions (see 
Pi), we obtain that a ~ (t^ — t)~ , i.e. we find a universe with a final Big Rip singularity. Thus 

~ (ts — ^ 00 at finite future time tg. 
Coming back to the case under consideration governed by eqn. ()3U() it follows that the effective 
value of w (for n / — 1) is 

1-n 2 



Weff 



n+l 



-1 + 



n + l 



(32) 



The reliability conditions one has to assume are respectively > and (n — 1) a > 0, which 

follow directly from (|30j) . Models for accelerating universe occur consequently for n>2orn<— 1. 
The first case corresponds to an effective quintessence when — 1 < Wgff < — | and one deals with 
an initial Big Bang type singularity at the origin. In contrast, the second case leads to an effective 
phantom {weff < — 1) (without its explicit introduction) with a Big Rip type final singularity (for 
a discussion of it in modified gravity see Q]), with the asymptotic behavior a{ts) ~ + 00. 
The case n = — 1 should be treated separately. In this case a = 0; this implies that the trivial so- 
lution R = = (j) emerges and we obtain a non-expanding universe model with a = oq = constant. 
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1. The radiating universe 



One can also combine the above case F{R) = R with the presence of radiation. Adding radiating 
matter {w = ^) gives rise to the following modified Priedmann equation: 

I I ^ +^ + ^ (33) 

\a 2hl Gbf{R)ae 3b 3ba^ ^ ^ 

To obtain exactly solvable models it is again convenient to specify to the case f{R) = aR^; we 
thus obtain {n ^ —1, ^,1): 

and we are able to compute the deceleration parameter to be 

3^ [(n - 1) g/g^] ^ g-^ + |fe^r?a-^ 

^ = 2iP 

A detailed analysis of this last equation is more complicated but the final answer is similar to 
the previous case: accelerated universe can be obtained only for n > 2 or n < —1.^ We stress, 
moreover, that the asymptotic value of the deceleration parameter for big enough radius a{t) in 
both cases remains the same, i.e. 

_ 2 - n 

Qasymp — ', T 
n + 1 

as in the radiation- free example (|28j) . Because of this the singularity types should be the same as 
well. For n < —1 (which, by the way, corresponds to the case 7 < 0), besides the final Big Rip 
singularity one can encounter the initial Big Bang singularity caused by radiation. To be more 
precise this follows from the qualitative behavior of (see (|34|) ): 

~ Aa-^^ + Ba-"^ (n < -1) 

where A and B can be simply deduced from (|34|) . In the case 7 < the first term dominates in late 
time universe, giving rise to a Big Rip like singularity. On the contrary, in the same case 7 < 0, we 
get that small values of a correspond to a behavior of ~ Ba'"^, which implies the presence of a 
Big Bang singularity at early time universe. In the case n > 2 (7 > 0) one obtains an accelerated 
and expanding model, such that the initial Big Bang is still present. Finally we remark that for 
n > 2 there exists a critical value for the radius 

n+l 

"2(n-l) " 



''crit 



n 



[(n-l)a/32] (36) 



^ Notice that some cases are excluded by the condition > and others give rise to deceleration. 
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which provides a transition from deceleration to acceleration. 

To summarize we have thus proven that, also in the Palatini formalism, non- linear gravity- matter 
systems based on the Lagrangian with F{R) = R lead to accelerating universe in the same way 
as it happens in the metric formalism. Big Rip singularities appear in the case of effective phantom 
models in vacuum universe, while quintessence models contain Big Bang like singularities. We have 
moreover shown that in the case of radiation models, where a radiation-like fluid is additionally 
considered, we have both Big Rip and Big Bang singularities appearing. 



A full set of possible FRW cosmological solutions of the above non-linear gravity-matter 
system based on the Lagrangians with F{R) = R has been found, under the specific hypothe- 
ses f{R) ~ i?", V{(j)) = and K = both in the case of vacuum and radiating universes. It turns 
out that the same solutions were already found and studied, for the same class of Lagrangians, in 
the framework of the purely metric formalism P| where they however represent a proper subset 
in the set of all possible solutions. We think that a bit of discussion is here appropriate on these 
arguments. 

It is well known that, in the matter free case L^^ = 0, L^at = 0, i.e. for pure -F(i?)-gravity, 

solutions of the Palatini formalism represent in g eneral a subset in the set of all possible solutions 

allowed by the purely metric formalism; see As a matter of fact in that case the conformal 

factor b = F'{R) relating the two metrics and becomes constant by virtue of the master 

equation; this implies that the scalar R is equal to the Ricci scalar of the metric g. As far as 

the field equations in the metric formalism are concerned, it implies also that the second order 

fl 

differential operator V^Vi, — g^uO acting on F'{R) vanish identically; see e.g. This finally 

implies in turn that the field equations in the metric formalism, expressed in the metric g^u, 
almost reproduce the field equations of the Palatini formalism. 

It is easy to check that exactly the same mechanism acts in our case, provided that f{R) ~ i?". As 
a consequence we obtain that, in the particular example we have here analyzed (i.e. F{R) = R), 
the Palatini formalism provides at the same time solutions to the metric one. 



B. First-order scalar F {R)-f (i?) non-minimal coupling 

We generalize hereafter the analysis to the more general family of models already introduced in 
(|H) for non-linear gravity-scalar system. To obtain exactly solvable models, which are preferable as 
they allow us to discuss relevant physical consequences, we will specify step by step the necessary 
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conditions ( / (R) = aR^ (a / 0; n / 1), r = and K = 0). We thus consider the more general 
class of Lagrangians: 



L = ^{F{R)+f{R) Ld)+Lmat 



(37) 



where = —^gf^'^d^4'd,y(j) + V {(/)) and Lmat represents a matter Lagrangian but now F{R) ^ R. 
We get now that the r.h.s. of the generalized Friedmann equation (jlSf) becomes: 



a b 

+ 



+ 



K 



c ("iw + l)r? 

+ TT7 + 



\a 2b a2 66 /(i?)a6 36 66a3(i+'") 
We proceed by specifying the model to the concrete particular case (n 7^ 1, m 7^ 2) 

/ (i?) = aR^ 
F [R) = LoR"" 



(38) 



(39) 



without matter (i.e. r = and tJ^ = 0) which leads again to an exactly solvable class of models. 
Combining the master equation with the on-shell value f{R)'^Lii = 0^a~^ (see ()17() ) it is simply 
possible to obtain the Ricci scalar as a function of the radius a: 



where the coefficient ^ is given by ^ 



R = 

. In this case it follows that: 



(40) 



Lu(2—m) 



1 

n-\-m 



6 = (mwC"""^ + nQ-i/32^-("+i)) a~^^^ 



6(m-l) 

60a 



6m 
m+n 



(41) 



Co a 



Since now 6 = 60a ^ is proportional to a power of a, it follows that the generalized Hubble 



parameter = ^ + ^ is, up to the constant factor 1 — '^^^j^^' , proportional to H. Therefore one 
has (again in the solvable case K = 0): 



3(m-l) 



ij2 ^ 



'66oa^" 36o 



Co 1 6_ 6_ 



(42) 



1__3_ 

This further implies that a ~ a ™+" (which resembles the case m = 1, already treated in Section 
Inj) . Consequently: 



a{t) ~ (t, ±t) — 



The deceleration parameter turns out to be q 



-1 + 



As for the very particular case, 



already considered in Section HTll we obtain by a similar analysis an expanding universe in the case 
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m + n > 3, presenting a Big Bang like singularity; conversely, for m + n < the expanding solution 
leads to a Big Rip finite time singularity. 

We stress that now the Palatini formalism and the metric one are completely different both from 
a qualitative and from a quantitative viewpoint. This is related to the fact that in this case the 
conformal factor b as given by © is no longer constant. Moreover we have seen that in the 
Palatini formalism field equations are easily solvable also in this more general case, while in the 
metric formalism this is by no means true. 

1. The case F{R) ^ R + fiG{R) 

We can consider the interesting case when F{R) = R + i^iG{R), which reproduces the case already 
treated in subsection IIII Al in the case /x = (or in the case G{R) = 0). The corresponding 
Lagrangian is then: 

L = ^{R + fiG{R) + f{R)Ld)+Lmat{'^) (43) 

and the master equation is: 

R + 2/iG (R) - fiG' {R)R + T= (/' {R)R-f (R)) L^-f (R) V (</<) (44) 

We specialize moreover to the case f {R) = ai?", G (R) = R^^, V{(l)) = and for simplicity 
Lmat = (the case of vacuum universe). Equation© gives in this case: 

b=l + f'{R)L, + ^^G'{R) = l + ^,=^ 

(45) 

c=^{R + f{R)La + fiG{R)) = Rj^^ - ^^lG (R) 

where the value of the conformal factor b is again constant, as much as /' (R) + f^G' (R) = 6—1 
is constant too. This implies that the solutions of the Palatini formalism are at the same time 
solutions of the metric formalism, as much as in the case F{R) = R previously examined (see our 
discussion after formula (|36|) '). It is however clear that in this case the Palatini formalism is much 
easier to handle and more effective then the purely metric formalism (even if it provides the same 
solutions). Substituting the actual value of along solutions into the master equation, as given 
by H17() . we obtain now: 

/?"+^ + 2^(1 - n)i?3« = a (n - 1) f3^a-^ (46) 

Since R is positive, in order to get an approximate solution of the polynomial equation above, one 
needs to distinguish two cases: 
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1. Small values of R {R « 1) and n >> or R >> 1 and n « lead to 

6 

R a "+i 

2. Big values oi R» 1 and n>>Oori?<<l and n << lead instead to 

i?~a~" (47) 

In any case we obtain a power law R ^ a^^, where A = in the first case and A = ^ in the 
second one. 

Owing to the fact that b is constant it is now easy to obtain the modified Hubble equations from 
H38() (in that case g is obviously chosen to be a FRW metric for cosmological applications): 

6(2n-l) 6(2n-l) 6 (2n - 1) ^ ^ 

which further re- writes as: 

with suitable chosen constants A, B,C, D. In the first case, recalling that A = — tt, one has: 

H'^ Ka-'^ + {A + B)a~'^ + Ca~^ 

6 „_ 

For late universe (a >> 1) the term a "+i dominates and the deceleration parameter is q = 
This gives the same qualitative behavior already found in the case F{R) = R previously discussed. 

12n 

For early universe (a << 1), instead, the term a "+i dominates, providing a Big Bang initial 
singularity. 

Correspondingly, for R a « the Friedmann equation reduces to 

H'^ ~ -Ka-^ + Ba-n +{A + C7)a-^ 

and consequently the term a « dominates in late universe. However this contradicts to our initial 
assumptions i? >> 1 and n>>Oori?<<l and n « for that case; see eqn. (|47|). For early 
universe, conversely, the term dominates and we obtain a universe fulfilled by radiation with 
a Big Bang initial singularity. 

IV. THE COSMOLOGICAL CONSTANT PROBLEM 

One of the alternative proposals to explain the present acceleration of the universe is related 
to dynamical cosmological constant models. This approach, however, has a lot of unsolved issues 
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related to the fact that there is no known mechanism, up to now, that guarantees zero or nearly zero 
energy in a stable energy configuration for the universe. Using cosmological constant-like models, 
we have to face with some problems concerning: (i) the small amount of energy of the vacuum, 
which is much smaller than we estimate it to be (the so-called cosmological constant problem); 
(ii) the nature of the dark energy which seems to dominate the universe and; (iii) the coincidence 
rohlem between the actual density of dark energy in the universe and the actual matter density 



It has been recently proposed in 3] that the apparent value of the cosmological constant can be 
alternatively determined by dynamical considerations in the framework of alternative theories of 
gravity. The true value of the vacuum energy is not zero, but settles down to a nearly zero value 
when the curvature of the universe reaches nearly zero values. This mechanism is produced by 
means of a feedback process in which the scalar equation depends on curvature, in such a way that 
the divergent coefficient of the kinetic term of the scalar field forbids the lowest energy state to be 
reached. The potential of the scalar field is chosen to have a minimum at a generic negative value 
for the vacuum energy. However, the singular kinetic term of the scalar field makes the dynamics 
of the scalar field to stop at zero, where the field itself gets frozen. (Of course, some fine-tuning of 
dynamical Lagrangian is assumed). 

This theory has been proven to have some very important properties, i.e. it is stable under radiative 
corrections and it leads to stable dynamics. No fine-tuning of the potential is necessary in this 
case, neither the introduction of an anthropic principle. 

The phenomenological theory introduced in ^ slightly differs from the one we studied in (pQ). The 
Lagrangian density is: 

L = ^(^ + aR^ + f{R)Lk^n-V{(p)] (50) 



where Lfcm = ~\9^^dii'P is the kinetic term of a scalar field Lagrangian and f(R) is a function 
of the Ricci scalar, which is assumed to be divergent at = 0. It should also be noticed that (|5()|1 
coincides with (HJ) for q = 0, F{R) = R and V{4>) = 0. 

We want here to examine if the same mechanism, generating a dynamical cosmological constant 
with suitable properties, works in the first order formalism. Let ghe a. flat FRW metric ©. First 
of all, even in the Palatini formalism the same_field equations hold for the scalar, i.e. they have 
the same form as ones in the metric formalism 

I [H{R)\ + ^H^fiR) + V'{ct>) = (51) 
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As far as field equations for the metric field g are concerned, we obtain in analogy with © that: 

(2^ + Vac^R^ (r,.u - + /' (R) LkinR^u - \f (R) Lking^. - \v{(p)g^,u = (52) 

and the scalar- valued structural field equation, controlling the solutions of (|52j) . can be simply 
obtained by taking the trace of (|52j) itself. This implies that, considering (p to evolve slowly, we 
have that T^((/>) can be approximated as a linear function, see P]: 

V{(fyc^cK-^{<l,-ct)o) (53) 

and the asymptotic behavior of cj) is consequently: 

K^^f{R) ~ -CK-^H-^ (54) 

This expression follows from the consideration that the term is nearly zero at small energies. 
From field equations H52() . in the approximation that the kinetic term is small compared to the 
potential term and assuming ~ 0, one obtains in analogy with Q|: 

m"^ ~ (55) 



This is exactly the same result already obtained in |3] in the metric formalism, proving that field 
equations from the Lagrangian H5U|) have the same limit at low energies, both in the Palatini 
and the metric formalism. It follows (see [3] for details), that at low energies (j) stalls at nearly 
y = 0, where V includes all the contributions to the cosmological constant. This consequently 
implies that the value of the cosmological constant stalls at a very small but non-zero value, in 
the limit of small curvatures of the metric 5, i.e. at small energies for the gravitational field. It 
has also been proven in P| that this behavior is stable under the effect of radiative corrections and 
leads to stable dynamics, provided that the minimum of the potential always occurs at a negative 
value. Hence, the same dynamical mechanism to solve the cosmological constant problem works in 
Palatini formulation. 



V. CONCLUSIONS 



In summary, the Palatini formulation of non-linear gravity-matter system has been developed. 
Using a scalar field Lagrangian as matter it is shown that such a theory provides the effective 
quintessence or effective phantom cosmology at late times in the same qualitative way as in the 
metric version. In addition, the gravity assisted dark energy dominance mechanism also occurs. 
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It is shown that an account of radiation maybe done with possible emergence of accelerating 
cosmology again. It is interesting to notice that in this case the emerging cosmology may contain 
Big Bang at early times and Big Rip at late times (for general review of late-time singularities, see 
ref and references therein). Eventually, the account of quantum effects near to singularities is 
necessary to see if the singularities are realistic or not. It would be really interesting to perform 
such a study within the Palatini approach. 

Any Palatini model of gravity possess its purely metric counterpart which leads to fourth order 
field equations. An advantage of Palatini formulation rely on second order field equations which 
turns out to be more easy to solve. In some very specific situations both approaches give rise to 
the same solutions but in general it is not true. The same occured here in the case of cosmological 
application. Although, one can expect that metric formulation may also produce accelerated solu- 
tions thus, however, it is much more difficult (if possible at all) to find them explicitly out. In this 
sense the Palatini formalism is more easy to handle and simpler to analyse then the corresponding 
metric formalism. 

It is obvious that any resonable model of gravity should satisfy the standard solar system tests. 
In the context of modify gravity this problem has been studied recently by several authors [2^ (both 
in the metric and Palatini approaches). It has been shown that, in principle. Palatini formalism 
provides a good Newtomian approximation. Particularly, the Schwartzschild limit of the theory 
can be recovered for positive powers n of the scalar R, greater than one n > 1. In the case of 
negative powers one should expect de-Sitter or anti-de-Sitter black holes instead [21 1. 

It has been also explicitly demostrated that dynamical mechanism to resolve the cosmological 
constant problem suggested in metric formulation works also in Palatini formulation (within the 
same choice of scalar-tensor Lagrangian). Again, it is very interesting to understand the role of 
quantum effects (as this is exactly the region where they are expected to be essential) in the above 
mechanism. However, to do so one needs to evaluate quantum effects in the Palatini formulation, 
something that is rather far from being a trivial task. 

Finally, in this paper we have limited ourselves to the consideration of non-linear scalar-gravity 
theories in the Palatini formulation. However, it is clear that non-linear spinor-gravity system 
(where the spinor part can be constructed from NJL-like model coupled to gravity ^22]) maybe 
more interesting from a cosmological point of view. Indeed, unlike the covariant derivative of a 
scalar which is nothing but a partial derivative, the covariant derivative of a spinor includes instead 
the connection, what definitely leads to new features in the spinor analog of Eq.7. This will be 
discussed elsewhere. 
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APPENDIX A: SOME USEFUL RELATIONS IN THE PALATINI FORMALISM 

We include hereafter some useful relation in the Palatini formalism for the benefit of the general 
reader. The Palatini formalism can be considered as a metric-affine theory for relativistic theories. 
Metric and connection are firstly considered as independent objects. We remark however that, 

without needing any reference to a metric, a connection can be introduced from which one can 
construct covariant derivatives. In order to produce covariant derivatives which transform as true 
tensors, under a coordinate transformation the connection itself needs to transform as: 

F'' + 



dxP dx'i^ dx'" dxP dx'P'dx"' 

Then for such a connection the quantity 

7->A fiK p)j -p^ -pjj -pX 

will then be a true tensor, regardless of whether the connection is or is not related to the metric 
in the standard metric based Christoffel symbol way 

However, noting that the metric based T^^ (g) also transforms as 

dx'^ dx'^ dx" . . dx'^ d'^xP 



(q) — _ — _ — _pp / N _|_ 



M^^v^/ dxPdx'P-dx'" dxP dx'P-dx'" ' 

we see that the difference between the general connection and the Christoffel symbol connection 
transforms as 



(q) ^ ^^^\rp -vp (q)] 
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with the difference thus being a true tensor. Then, with covariant derivatives being Unear in 
the connection, the difference between the quantity T-^^ (where T^'^ is the stress-energy tensor of 
matter), as calculated with Tra, and T-^ , as calculated with Traig), is itself a true tensor, and 
since T-^ as calculated with Tra is necessarily a true tensor, it follows that T-^ as calculated with 
Trcr{g) is one as well. It is consequently meaningful to take the covariant derivative of equation 
of motion ((2j with respect to the Christoffel symbol based connection: it transforms like a true 
tensor. 
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